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Electronic reconstruction at the polar interface LaMn03/SrMnC>3 (LMO/SMO) (100) resulting 
from the polar catastrophe is studied from a model Hamiltonian that includes the double and super 
exchange interactions, the Madelung potential, and the Jahn- Teller coupling terms relevant for the 
manganites. We show that the polar catastrophe, originating from the alternately charged LMO 
layers and neutral SMO layers, is quenched by the accumulation of an extra half electron per cell 
in the interface region as in the case of the LaA10;i/SrTi03 interface. In addition, the Mn e g 
electrons leak out from the LMO side to the SMO side, the extent of the leakage being controlled 
by the interfacial potential barrier and the substrate induced epitaxial strain. The leaked electrons 
mediate a Zener double exchange, making the layers adjacent to the interface ferromagnetic, while 
the two bulk materials away from the interface retain their original type A or G antiferromagnetic 
structures. A half-metallic conduction band results at the interface, sandwiched by the two insulating 
bulks. We have also studied how the electron leakage and consequently the magnetic ordering are 
affected by the substrate induced epitaxial strain. Comparisons are made with the results of the 
density-functional calculations for the (LMO)e/(SMO)4 superlattice. 

PACS numbers: 75.70.Cn; 73.20.-r; 75.70.-i 
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I. INTRODUCTION 

Polar interfaces have commanded considerable inter- 
est recently because extra electrons may migrate to the 
interface to "heal" the polar catastrophe and these in- 
terfacial electrons may exhibit unusual two-dimensional 
properties. A polar interface of current interest is 
LaA10 3 /SrTi0 3 (LAO / STO)i^££ , which shows a 
variety of phenomena such as the Kondo resistance min- 
imum, superconductivity, magnetism, and metallic or in- 
sulating behavior under varying circumstances* 2 - 

This paper is devoted to a theoretical study of the in- 
terface between LaMn03 and SrMnOs, which is a polar 
interface (Fig. [T]) with the extra twist that both con- 
stituent materials are also magnetic. Presence of the Mn 
localized moments introduces a new interaction channel 
for the electron gas forming at the interface. In fact, 
in a similar system, viz., a single 5-doped LMO layer 
in a SMO matrix, where an extra electron per cell be- 
comes introduced at the £-doped layer, the formation of 
a spin-polarized two-dimensional electron gas was pre- 
dicted due to the Zener double exchange between the itin- 
erant electrons and the localized Mn moments. 7 On the 
experimental front, it has been demonstrated that high- 
quality superlattices of the manganites can be grown by 
MBE and the magnetic properties can be controlled by 
the superlattice period as well as by substrate-induced 
strain] 8 ' 9 : 10 : 11 : 12 : 13 : 14 Various experimental and theoreti- 
cal studies have found diverse magnetic phases for these 
superlattices^ i 1 ^a2 i i3a4 i i5 i 6 i 7 1 i8 i 9^o 1 2i 

Both LMO and SMO are magnetic with LMO being 
a type A antiferromagnetic (AFM) insulator with the 
Mn-ifge 1 configuration and consisting of alternating pos- 
itively charged LaO and negatively charged MnO layers 
along the (100) direction. On the other hand, SMO is 
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FIG. 1: Illustration of the polar catastrophe for the 
LMO/SMO (100) interface. The layers are neutral on the 
SMO side, but are charged ±1 on the LMO side leading to the 
unrestricted growth of the Coulomb potential (dashed line) 
and the scenario where it is healed by accumulating 0.5 elec- 
trons per cell on the interfacial Mn02 layer (solid line). The 
layer numbers indicated in this figure are used throughout the 
paper. 



a type G AFM insulator with Mn-ij^e configuration 
with neutral SrO and Mn02 layers. As a result of the 
layer charge configuration, a polar catastrophe - the di- 
vergence of the Coulomb potential away from the inter- 
face - arises in this system (Fig. [T]) as in the case of the 
LAO/STO interface. We note that unlike the latter case, 
where different layer termination produces two different 
n and p type interfaces, for the LMO/SMO interface, the 
Mn02 layers being common to both, we just have a sin- 
gle type of interface, which by simple electron counting, 
is expected to be n type, with extra electrons coming to 
the interface region as illustrated in Fig. [TJ We study 
the issues of polar catastrophe, charge leakage, and mag- 
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FIG. 2: Layer projected Mn e g occupancy across the interface, 
as obtained from the model Hamiltonian with dielectric con- 
stant e = 10. About half an electron per Mn atom accumulate 
at the interface to remove the polar catastrophe and, further- 
more, a small number of electrons leak from the LMO to the 
SMO side. Electrons at the interface form a spin-polarized 
2DEG, sandwiched between the two insulating bulks. 



netism at the interface by using a tight-binding model 
Hamiltonian that includes all the relevant interactions in 
the system, complementing the study with the ab initio 
density-functional calculations. 

The main results that emerge from our work are the 
following: (i) Half an electron per unit cell accumulates at 
the interface to avoid the polar catastrophe much like the 
case of LAO/STO, although now the electrons are much 
more confined to the interface (Fig. [2]). (ii) In addition, 
the Mn (e g ) electrons leak from the LMO to the SMO 
side, spreading to several layers, (iii) While the LMO and 
SMO layers away from the interface retain, respectively, 
the type A and type G AFM structures of the bulk, the 
interface becomes ferromagnetic (FM) due to the double 
exchange between the itinerant carriers and the Mn core 
spins, (iv) Finally, we show that the magnetic ordering 
at the interface is sensitive to the substrate-induced epi- 
taxial strain. However, our results suggest the absence 
of any canted magnetic state at the interface. 



II. MODEL HAMILTONIAN 

Density-functional calculations on the manganitcs 
have provided important insights into the electronic 
structure ! 22 ' 23 The important electrons near the Fermi 
energy are the itinerant Mn (e g ) electrons, which inter- 
act with the localized Mn (t2 S ) core spins via the double 
exchange . 24 ' 25 The other important terms to consider are 
the Jahn- Teller coupling of the e g electrons with the oc- 
tahedral vibrational modes and the Madelung potential 
of the constituent atoms that would lead to the polar 
catastrophe without the electronic reconstruction at the 
interface. 

Incorporating these key interactions, we construct the 



model Hamiltonian 

H = y^\e ia + - ^ Mi V g v )c\ a c ia 

ia V 

+ t ^,3pcos{e ij /2){c\ a c jf} + H.c.) 

{ij)a0 

+ f E™ 2 + E^ + ^E^-^' w 

i i (ij) 

which describes the motion of the Mn (e g ) electrons in a 
matrix of Mn core spins. These electrons are effectively 
spinlcss as discussed later. The e g electrons are restricted 
to the Mn sites, but all atoms in the structure (Mn, Sr, 
La, and O) contribute to the electrostatic potentials that 
these electrons see, so that the Madelung matrix My goes 
over all atoms. In the Hamiltonian above, cj creates an 
electron at the i th Mn site and a is the orbital index of the 
e g electron (x 2 — y 2 or z 2 — 1), rii a is the corresponding 
number operator, (ij) denotes summation over nearest 
neighbor pairs, and while the i summation runs over only 
the Mn sites, the v summation runs over all atoms in the 
structure. 

The first term in the Hamiltonian describes the onsite 
energy, where the energies Si a of the two e g orbitals could 
be split due to strain (considered later in the paper) . The 
Madelung energy is given by E Ma d = ^ Yl, V u' M wq_ v q_v', 
where q„ denotes the total charge (ionic + electronic) 
of the i^-th atom, e is the dielectric constant, the di- 
agonal terms of the Madelung matrix M vu t exclude the 
Coulomb contribution from the same site, and the fac- 
tor of two comes from double counting. The Madelung 
potential seen by the e ff electron is given by Vuad — 
dEMad/driia = \ ^2, v Mi V q Vl which appears in the first 
term in the Hamiltonian. The Madelung matrix for the 
structure is obtained by using the standard Ewald sum- 
mation method. 

The second term in the Hamiltonian is the electronic 
hopping energy. We have taken the Hund's energy Jjj to 
be oo, so that the coupling between the core and the itin- 
erant spins Ti-Hund = — Jh J2ia Si'Sia makes the electron 
state inaccessible, where the electron spin is anti-aligned 
with the local t^g core spin. Thus the itinerant electrons 
are effectively spinless and the electron hopping is di- 
minished by the Anderson-Hasegawa cos(#y/2) factor—, 
where By is the angle between the (classical) core spins 
at the two sites. The nearest-neighbor hopping integral t 
depends on the relative positions of the two Mn sites in 
the lattice. For nearest-neighbor hopping in the xy-plane 
or along the z-axis, we have: 

e v-Y^( 1 -V3\ / 1 0\ () 

*a/9 - 4 ^ 3 J - V-"^O0j' [ > 

where V a is the dda matrix element and we have ne- 
glected the much smaller dd5 interaction. The third term 
in the Hamiltonian is the Coulomb repulsion between the 
two e g orbitals on the same site. 
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The Jahn- Teller (JT) coupling term on each site is 
given by 

Hj T = \k{Q\ + Ql) - 9 (Q 3 t z + Q 2 t x ), (3) 

where f is the pseudospin describing the two e g orbitals, 
viz., | T) = \z 2 - 1) and | |) = \x 2 - y 2 ). With the 
corresponding creation operators being cjj and c\ 2 , re- 
spectively, we have: r = J2 a 8 c \oJu(}Cif3, where the greek 
indices denote the orbitals and i is the site index. The 
quantities Q2 and Q3 are, respectively, the basal plane 
distortion mode and the octahedral stretching mode at 
the i-th site, K is the elastic stiffness constant, and g is 
the JT coupling strength. The final term in the Harnil- 
tonian describes the AFM superexchange interaction be- 
tween the Mn t2 fl core spins. 

The typical values of the Hamiltonian parameters used 
in our calculations are: K = 9 eV/A 2 following from the 
optical studies on Lao.ssSro.isMnO^, g = 2.0 eV/A, 
and V<j = -0.5 eV following earlier density- functional 
results 2 -, J = 26 mcV, estimated from the Need temper- 
ature of CaMn0 3) 28 i 29 a compound similar to SrTiC>3, 
and finally, U = 3 eV, which is a reasonable for the 
Mn (e s ) electrons. The values for the JT distortions in 
the LMO bulk are^ Q° = 0.28 A and Q° = -0.10 A. 
These distortions are bulk like on the LMO side and go 
to zero on the SMO side, where there are no e g elec- 
trons. There is a transition region in between, where, 
for simplicity, we have taken the distortion strengths to 
be linearly dependent on the site electron occupancy of 
the e g orbitals, which would be the result for the iso- 
lated octahedron. This is easily seen from an examina- 
tion of the single-site JT Hamiltonian Eq. [3] If the elec- 
tron occupancy n is taken to be a continuous variable 
as appropriate for a mean-field model, it immediately 
follows from the diagonalization of the Hamiltonian ([3]) 
that the optimized magnitude of the distortion is given 
by: (Q2 + Q3) 1 / 2 = gn/K. Following this argument, 
we use the distortion magnitudes at each Mn site as: 
Q2 = n x Q^, where is the magnitude for bulk LMO 
and similarly for Q3. 

There is one more point to be made regarding the 
model. The measured value of the dielectric constant for 
LMO is 16-18, while for SMO, it is roughly 35. 31 ' 32 From 
work on the LAO/STO interface, it has been argued that 
the dielectric constant for STO is drastically reduced 3 - 3 - 
in the presence of an electric field, which exists at the in- 
terface. We however use a uniform dielectric constant in 
our model calculations following earlier authors,— ^ which 
is reasonable within the spirit of our model. We find that 
the electron density profile near the interface region ob- 
tained from our model using e 1=3 8 — 10 fits well with 
the ab initio density-functional results. Unless otherwise 
stated, we have taken e = 10 in our model calculations. 

We solve the model Hamiltonian within the Hartree- 
Fock mean field approximation. A supercell geometry is 
adopted for the convenience of calculation; we find that 
the (LMO) 10 /(SMO) 10 cell is large enough for our pur- 



TABLE I: Relative energies (per Mn interface atom) for dif- 
ferent magnetic configurations. The intra-layer magnetic cou- 
plings are denoted by the symbols F or A (FM or AFM or- 
dering within the plane), while bracketed symbols (F) or (A) 
denote the coupling between the adjacent interface layers. (F, 
A) means that half the interlayer bonds are FM and the other 
half are AFM. The lowest energy structure, corresponding to 
the first line of the Table, is shown in Fig. [3j 
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poses. For a fixed core spin configuration, beginning with 
an initial set of Mn site occupancies, the band structure 
was computed in the Brillouin zone, from which new site 
occupancies were determined, and the process was re- 
peated until self-consistency was reached. 

To compare with the results of our model calculations, 
we have also performed an ab initio density functional 
calculation using the linear muffin-tin orbitals method 
(LMTO) with the generalized gradient approximation 35 
to the exchange-correlation functional with the Coulomb 
correction (GGA + U) using a somewhat smaller super- 
cell, viz., (LMO) 6 /(SMO) 4 , and a value of U = 3 eV, 
which is reasonable for the Mn (e g ) electrons. 



III. RESULTS AND DISCUSSIONS 

We examined the computed total energy for different 
configurations of the core Mn spins and found that the 
minimum energy structure is one where away from the 
interface, the G and A type antiferromagnetism of the 
bulks are recovered, while the interface region becomes 
ferromagnetic. The magnetic configurations at the in- 
terface are characterized by magnetic ordering within 
each plane, indicated by unbracketed symbols in Table 
I, and between the adjacent planes, which are indicated 
by bracketed symbols. Deeper in the bulks, both LMO 
and SMO are AFM type A and G, respectively. The 
computed energies for different magnetic configurations 
are listed in Table I and the configuration with the low- 
est energy is illustrated in Fig. [3] The results presented 
in the following sections correspond to the lowest-energy 
configuration, unless otherwise stated, and we note that 
quantities such as charge reconstruction and Madelung 
potentials are rather insensitive to the magnetic configu- 
ration. 



FIG. 3: (Color online) Minimum-energy magnetic configu- 
ration at the LMO/SMO interface for typical parameters 
and unstrained condition as obtained from both the density- 
functional and model calculations. The first MnCh layer on 
the SMO side can be tuned FM or AFM by changing the 
strain condition as discussed in the text. 
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FIG. 4: Madelung potential seen by the Mn ions in various 
layers without and with the charge reconstruction as obtained 
from the self-consistent calculation using the model Hamilto- 
nian (Eq. 1). Without the charge reconstruction, the poten- 
tial grows unrestricted away from the interface (polar catas- 
trophe), which is healed by accumulation of half an electron 
per cell at the interface. The reconstructed charge at the in- 
terface has a dipole moment that leads to a potential step 
(bottom figure). 

A. Polar catastrophe and Interfacial charge 
reconstruction 

As widely discussed in the context of the LAO/STO 
interface, alternately stacked positive and negative lay- 
ers - in the present case LaO and MnC>2 layers respec- 
tively - lead to a divergent Coulomb potential, the so 
called polar catastrophe. In Fig. [3J we have plotted the 
initial Madelung potential of the charged layers (taking 
the nominal charged states for all atoms) and the final 
potential after charge reconstruction has occurred and 
self-consistency has been achieved. 



-3-2-10 1 2 3 
Layer No. 

FIG. 5: (Color online) Mn e 9 occupancy across the 
LMO/SMO interface for different values of the dielectric 
constant e (symbols connected by lines) compared with the 
density-functional results for the (LMO)e/(SMO)4 superlat- 
tice (filled squares). An increased e results in a diminished 
Madelung potential and, consequently, electrons leak deeper 
into the bulk. 



Charge reconstruction at the interface occurs in two 
ways as may be seen from Fig. [5] First, a monopole 
charge, half electron per cell, or 3.4 x 10 14 e/cm 2 accumu- 
lates at the interface region to cancel the original electric 
field discontinuity at the interface. A significant portion 
of these electrons reside on the interfacial MnC>2 layer as 
seen from the figure. The interface region demands this 
monopole charge in order to quench the polar catastro- 
phe. In our supercell geometry, which has two identical 
interfaces per supercell and also an extra electron per cell 
after satisfying the nominal ionic charges of all atoms, 
this electron was quite conveniently shared between the 
two interfaces. Thus we have half an electron per cell 
available at each interface, exactly the amount needed to 
quench the polar catastrophe. In the experimental situa- 
tion, this monopole charge must come from the surfaces 
and/or bulk defect states as necessary, so that the de- 
mand for the monopole charge at the interface can be 
met. A second aspect of the charge reconstruction is 
that the Mn (e fl ) electrons leak from the LMO side into 
the SMO side, creating a dipole moment at the interface, 
which results in a potential discontinuity of about 2 eV 
as seen from the bottom part of Fig. 2] The magnitude 
of this surface dipole moment density is estimated to be 
p m O.lle/A. 

The interfacial electronic reconstruction is consistent 
with the layer projected densities of states (DOS) shown 
in Fig. [5] Here we see that at the interfacial M11O2 
layer, approximately one fourth of the e g states are occu- 
pied indicating a net occupancy of about half an electron 
per Mn atom. The DOS for the layers away from the in- 
terface resemble the bulk electronic configuration of the 
respective materials. Layers close to the interface have 
either a small number of electrons in the e g bands (SMO 
side) or a small number of electrons missing (LMO side), 
consistent with the electron distribution of Fig. [5] 
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FIG. 6: (Color online) Layer projected Mn e g densities of 
states obtained from the model Hamiltonian indicating extra 
electrons or holes accumulated at the interfacial layers, while 
away from the interface, the e s states are either half filled 
(LMO) or unoccupied (SMO). The electrons are effectively 
spinless, with only the spin state aligned parallel to the Mn 
core spin at a particular Mn site allowed to be occupied (Jh ~ 
oo), so that we have a single spin channel in this figure. 



B. Density- functional results and half- metallic 
behavior 



In this subsection, we report the results of our density- 
functional theory (DFT) calculations, which validate the 
results of our model and also address the issue of half 
metallicity of the conduction bands at the interface. Half- 
metallic behavior, where one spin band is metallic, while 
the other is insulating, was already predicted for the delta 
doped LMO/SMO structure, where a single LMO layer 
is doped in a matrix of SMO.— There, the extra electrons 
coming from the LMO layer became confined in the elec- 
trostatic potential well of the LMO layer producing a 
ferromagnetic alignment of the core Mn spins at the in- 
terface and in turn became spin polarized owing to the 
Zeeman field of the Mn moments. Our DFT studies show 
that we have a similar situation for the LMO/SMO inter- 
face, with the difference that now the interfacial electrons 
originate from the polar catastrophe rather than from the 
La dopants in the delta doped structure. 




FIG. 7: Variation of the energy of the lowest Mn e 3 state of 
each Mn02 layer, obtained using DFT-LMTO from the layer- 
projected wave-function characters for the (LMO)e/(SMO)4 
superlattice. This variation in energy is indicative of the po- 
tential experienced by the Mn e 9 electrons across the inter- 
face. 



The DFT calculations were performed for the 
(LMO)e/(SMO)4 superlattice using the linear muffin-tin 
orbital (LMTO) with the gradient approximation for the 
exchange correlation functional and including the on-site 
Coulomb term (GGA+U) with U = 3 eV. Both the in- 
plane and out-of-plane lattice parameters were taken to 
be the average lattice parameters of LMO and SMO. 
Lattice relaxation effects were not included, which are 
expected to cause the electrons to spread somewhat fur- 
ther into the bulk away from the interface, but otherwise 
not change the essential physics of the problem. Effect of 
lattice relaxation for a similar system, viz., delta-doped 
SrTi03/(LaTi03)i/SrTi03 was discussed in our earlier 
work using detail density functional calculations.— 

In Fig. [7| we have shown the relative potential seen 
by the Mn e ff electrons at each Mn02 layer as obtained 
from the DFT calculations. The variation in the poten- 
tial across the interface is obtained by calculating the 
energy of the lowest Mn e g band state in each M11O2 
layer, which can be obtained from the layer-projected 
wave-function characters. From the figure, we see that 
the magnitude of the potential discontinuity at the inter- 
face (about 2 eV) is similar to the one obtained from the 
model calculation (Fig. 0|. 

The layer projected densities of states for the 
LMO) 6 /(SMO) 4 superlattice, as obtained from the DFT 
calculations, are shown in Figs. [5] and [5] From the fig- 
ures, we see that the DOS for the innermost Mn02 layers 
show bulk-like behavior, while the occupation of the Mn 
e ff states for the interfacial layers are very similar to the 
results obtained from our model calculation (Fig. [5]). 

A notable feature of the density functional results is 
that, for each layer the minority-spin states are com- 
pletely unoccupied at the Fermi level, leading to the 
half- metallic behavior. In other words the two dimen- 
sional electron gas formed at the interface is completely 
spin-polarized. We note that this behavior is obviously 
absent in the widely studied LAO/STO interface, which 
lacks any magnetic atoms. 
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FIG. 8: Density-functional results for the majority and mi- 
nority spin Mn-d densities of states for the SMO and LMO 
layers away from the interface indicating the bulk-like insu- 
lating behavior. Results are for the innermost Mn02 layers 
in the (LMO) 6 /(SMO) 4 superlattice. 
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in-plane 
Tensile 



FIG. 10: Energy splitting of the Mn e g orbitals due to in- 
plane compressive and tensile strain conditions used in our 
model Hamiltonian. 



Layer population in SMO 









: Layer #0 _ 






0.4 










o 












to 












Q. 
13 












O 

o 
o 


0.2 


^_J_ayer #1 






Layer # 1 
SMO 


<D 




Layer #2 











-2 -1 





1 2 



Compressive 



Strain (A) (eV) 



Tensile 



FIG. 11: Electron leakage to the individual Mn02 layers on 
the SMO side as a function of strain. 



C. Effect of strain on electron leakage and 
magnetic ordering 

In this subsection, we study the effect of strain on 
the electron leakage across the interface, which in turn 
affects the magnetism. Density-functional calculations 
have shown that strain alters the relative energy be- 
tween the two e s orbitals, causing a change in the orbital 
ordering^ A change of the symmetry of the occupied 
state from x 2 — y 2 to z 2 — 1 would increase electron hop- 
ping across the interface leading to an increase of electron 
leakage. Strain is taken into account in our model Hamil- 
tonian ((T|) via the on-site energy: 




A 



(a 



y 2 ) 
-D 



(4) 



where a positive A correspond to an in-plane tensile 
strain (see Fig. [T0|) . 

The computed layer occupancy of the electrons is plot- 
ted in Fig. QTJ where we have shown the electron leakage 
to various Mn layers in the SMO side. As expected, the 
layer occupancy diminishes with increasing A, a result 
of the lower on-site energy for the x 2 — y 2 orbital, which 
leads to a larger x 2 — y 2 character of the occupied e g elec- 
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FIG. 12: The total energy as a function of the canting angle 
8 and strain A. The parameter 8 is the angle between the 
neighboring Mn core spins in the first Mn02 layer on the 
SMO side. The figure shows that, the considered Mn02 plane 
either stabilizes with FM order or AFM order depending on 
the strain condition. No canted magnetism is found. 



trons and consequently to diminished electronic hopping 
across the interface. The layer occupancy of the interface 
layer # remains close to 0.5 e~, the density needed to 
quench the polar catastrophe. 

The magnetism in the manganites is determined by 
a competition between the superexchange between the 
core t2 g spins and the double exchange between the core 
spins and the itinerant e g carriers. The strength of the 
latter clearly depends on the concentration of the itiner- 
ant carriers. In SMO bulk, the AFM super exchange is 
the only term, since there are no e g electrons, which leads 
to a Neel type G order. As the concentration of the e g 



2.0 



1.0 



> 

< 



-1.0 



AFM 


Layer # 1 in SMO 


/ FM 





10 



20 



30 



40 



FIG. 13: Magnetic phase diagram for the first Mn02 layer on 
the SMO side as a function of the dielectric constant e and 
strain A. Rest of the layers maintain the magnetic behavior of 
the corresponding bulk materials, either AFM or FM within 
each layer, throughout the range of parameters considered. 



electrons is increased, the competition between the super 
and the double exchange could lead to a canted magnetic 
state, eventually resulting in a ferromagnetic state if the 
double exchange dominates^ 

To study if strain can affect the magnetic behavior by 
modifying the electron leakage, we have solved the Hamil- 
tonian with different orientation of the Mn core spins in 
order to obtain the ground-state magnetic configuration. 
Canted states were also considered in addition to the FM 
and AFM states. Within the range of strain studied, the 
strain being parametrized by the energy splitting param- 
eter A, we find that the magnetic ground state does not 
change from the one shown in Fig. [3] except for the mag- 
netic configuration of the layer # 1, which can be altered 
between FM and AFM depending on the strain condi- 
tion. In particular, the electron leakage to this layer can 
be large enough to produce a net ferromagnetic interac- 
tion. 

Fig. [12] shows the total energy as a function of the 
canting angle 9 between neighboring core spins in the 
MnC>2 layer # 1. There is a transition between the FM 
and AFM states, but we do not find the occurrence of a 
canted state within our model. Since the electron leakage 
beyond this layer is small, we see the bulk magnetic or- 
dering beyond the first layer. Fig. [13]shows the magnetic 
phase diagram for the same Mn02 layer as a function of 
strain and the dielectric constant. 



IV. SUMMARY 



In summary, by solving a model Hamiltonian in the 
Hartree-Fock mean-field theory, we studied the charge re- 
construction at the polar interface of LMO/SMO. The re- 
sults were complemented by ab initio density-functional 
calculations of the (LMO)g/(SMO)4 superlattice. Two 
types of electronic reconstructions were found at the in- 
terface. First, there is an accumulation of an extra half 
electron per cell in the interface region as in the case 
of the LAO/STO interface in order to quench the po- 
lar catastrophe. Secondly, the Mn e g electrons leak out 
from the LMO side to the SMO side and alter the mag- 
netism at the interface, while away from the interface, 
the magnetism of the respective bulk materials is pre- 
served. Our calculations suggest the presence of a half- 
metallic two-dimensional electron gas in the interfacial 
region under certain strain conditions. Indeed, experi- 
mental evidence for half-metallicity for this interface was 
recently obtained from magneto-optical Kerr effect ^1 We 
note that the LMO/SMO interface is a magnetic coun- 
terpart of the much studied polar interface between LAO 
and STO and as such presents an extra degree of freedom 
for the study of the two-dimensional electron physics. 

This work was supported by the U. S. Department of 
Energy through Grant No. DE-FG02-00ER45818. 
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